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Abstract 

In two remarkable recent papers, hep-th/0610248 and hep-th/0610251, the complete pla- 
nar perturbative expansion was proposed for the universal function of the coupling, f{g), 
appearing in the dimensions of high-spin operators of the AT = 4 SYM theory. We study 
numerically the integral equation derived in hep-th/0610251, which implements a resum- 
mation of the perturbative expansion, and find a smooth function that approaches the 
asymptotic form predicted by string theory. In fact, the two leading terms at strong 
coupling match with high accuracy the results obtained for the semiclassical folded string 
spinning in AdS^. This constitutes a remarkable confirmation of the AdS/CFT correspon- 
dence for high-spin operators, and equivalently for the cusp anomaly of a Wilson loop. 
We also make a numerical prediction for the third term in the strong coupling series. 
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1 Introduction 



The dimensions of high-spin operators are important observables in gauge theories. It is 
well-known that the anomalous dimension of a twist-2 operator grows logarithmically for 
large spin S, 

A-S^f{g)]nS + 0{S'), ^ = ■ (1) 

This was demonstrated early on at 1-loop order [1] and at two loops [2] where a cancellation 
of In^ S terms occurs. There are solid arguments that (1) holds to all orders in perturbation 
theory [3, 4], and that it also applies to high-spin operators of twist greater than two [5]. 
The function of coupling f{g) also measures the anomalous dimension of a cusp in a 
light-like Wilson loop, and is of definite physical interest in QCD. 

There has been significant interest in determining f[g) in the J\f = A SYM theory. 
This is partly due to the fact that the AdS/CFT correspondence [6] relates the large g 
behavior of / to the energy of a folded string spinning around the center of a weakly 
curved AdS^ space [7]. This gives the prediction that f{g) — > Ag at strong coupling. The 
same result was obtained from studying the cusp anomaly using string theory methods 
[8]. Furthermore, the semi-classical expansion for the spinning string energy predicts the 
following correction [9] : 

f(g)=Ag-^ + 0{l/g). (2) 

TT 

It is of obvious interest to confirm these explicit predictions of string theory using extrap- 
olation of the perturbative expansion for f{g) provided by the gauge theory. 

Explicit perturbative calculations are quite formidable, and until recently were available 
only up to 3- loop order [10, 11]: 

fi9)-Sg'-lnV + f^nV + Oig% (3) 

Kotikov, Lipatov, Onishchenko and Velizhanin [10] extracted the J\f — 4 answer from the 
QCD calculation of [12] using their proposed transcendentality principle stating that each 
expansion coefficient has terms of the same degree of transcendentality. 

Recently, the methods of integrability in AdS/CFT^ [16, 17], prompted in part by 
[18, 7], have led to dramatic progress in studying the weak coupling expansion. In the 
beautiful paper by Beisert, Eden and Staudacher [19], which followed closely the important 
earlier work in [20, 21], an integral equation that determines f{g) was proposed, yielding 
an expansion of f{g) to an arbitrary desired order. The expansion coefficients obey the 

^For earlier work on integrability in gauge theories, see [13, 14, 15]. 
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KLOV transcendentality principle. In an independent remarkable paper by Bern, Czakon, 
Dixon, Kosower and Smirnov [22], an explicit calculation led to a value of the 4- loop term. 



which agrees with the idea advanced in [22, 19] that the exact expansion of f{g) is related 
to that found in [20] simply by multiplying each ^-function of an odd argument by an i, 
C(2n+ 1) — > i({2n+ 1). The integral equation of [19] generates precisely this perturbative 
expansion for f{g). 

A crucial property of the integral equation proposed in [19] is that it is related through 
integrability to the "dressing phase" in the magnon S-matrix, whose general form was 
deduced in [23, 24]. In [19] a perturbative expansion of the phase was given, which starts 
at the 4-loop order, and at strong coupling coincides with the earlier results from string 

theory [23, 25, 26, 27, 21]. An important requirement of crossing symmetry [28] is satisfied 
by this phase, and it also satisfies the KLOV transcendentality priciple. Therefore, this 
phase is very likely to describe the exact magnon S-matrix at any coupling [19], which 
constitutes remarkable progress in the understanding of the A/" = 4 SYM theory, and of 
the AdS/CFT correspondence. 

The papers [19, 22] thoroughly studied the perturbative expansion of f{g) which follows 
from the integral equation. Although the expansion has a finite radius of convergence, as is 
customary in certain planar theories (see, for example, [29]) , it is expected to determine the 
function completely. Solving the integral equation of [19] is an efficient tool for attacking 
this problem. In this paper we solve the integral equation numerically at intermediate 
coupling, and show that f{g) is a smooth function that approaches the asymptotic form (2) 
predicted by string theory for g > 1. The two leading strong coupling terms match those 
in (2) with high accuracy. This constitutes a remarkable confirmation of the AdS/CFT 
correspondence for this non-supersymmetric observable. 

The structure of the paper is as follows. The integral equation of [19] is reviewed and 
solved numerically in section 2. An interpretation of these results and their implications 
for the AdS/CFT correspondence are given in section 3. 

2 Numerical study of the integral equation 

The cusp anomalous dimension f{g) can be written as [20, 19, 30] 




(4) 



f{g) = 16/a(0), 



(5) 



where a{t) obeys the integral equation 




(6) 
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with the kernel given by [19] 

K{t, t') = K^"^^ {t, t') + (t, t') . (7) 
The main scattering kernel K^"^^ of [20] is 

and the dressing kernel K^'^'' is defined as the convolution 

POO f" 

K^^\t, t') = 4/ / dt"K,{t, 2gt")^^Ko{2gt'\ t'), (9) 
Jo — 1 

where Kq and Ki denote the parts of the kernel that are even and odd, respectively, under 
change of sign of t and t': 

(+ +'\ tJi{t)Jo{t') - t'Jo{t)Ji{t') 2 ^ 
^Q[t,t) = ^2 _ ^ 2^[2n-l)J2n-l[t)J2n-l[t), (10) 

We find it useful to introduce the function 

s{t)^'^a{t), (12) 
in terms of which the integral equation becomes 

POO 

s{t) = K{2gt, 0) - Ag^ / dt'K{2gt, 2gt')^—s(t'). (13) 

Jo — 1 

Again, f{g) = 16^2^(0). 

Both and X'^^ can conveniently be expanded as sums of products of functions of 
t and functions of t': 

n oo 

K^"'\t,t') = Ko{t,t') + K^{t,t') = —Y,nUt)Ut'), (14) 

n=l 

and 

K^'Kt^t') = f: f: ^^^^^J^^2n,2„.-lJ2n(t)J2^-l(i')- (15) 
n=l m=l 

This suggests writing the solution in terms of linearly independent functions as 

s{t) = T.sJ-^. (16) 

n>i ^y'' 



so that the integral equation becomes a matrix equation for the coefficients s„. The 
desired function f{g) is now 

fig) = SgW (17) 
It is convenient to define the matrix as 

^„„.f d» "-'g>_y . (18) 

Using the representations (14) and (15) of the kernels and (16) for s{t), the integral 
equation above is now of the schematic form 

Sn^K-T. {k^^ + 2K(i) S^, (19) 



m>l 



whose solution is 



The matrices are 



^ h. (20) 



1 + /STM + 2Ki^) 



K^Zl = 2(7VZ)„^, (21) 
K^i = 2{CZ\^, (22) 
= 2{PNZQN)^^, (23) 

where Q = diag(l, 0, 1, 0, ...), P = diag(0, L 0. L ...), N = diag(l, 2, 3, ...) and the vector 
h can be written as h = (1 + 2C)e"'", where e = (1,0,0,...). The crucial point for the 
numerics to work is that the matrix elements of Z decay sufficiently fast with increasing 
m,n (they decay like e~™^("*'")/^). For intermediate g (say g < 20) we can work with 
moderate size d by d matrices, where d does not have to be much larger than g. The 
integrals in Znm can be obtained numerically without much effort and so we can solve for 
the Sn- We find that the results are stable with respect to increasing d. 

Even though at strong coupling all elements of Znm are of the same order in 1/g, 
those far from the upper left corner arc numerically small (the leading terms in 1/g are 
suppressed by a factor [(m — n — l){m — n + l){m + n — l){'m + n + for m ^ n±l). 
This last fact makes the numerics surprisingly convergent even at large g and moreover 
gives some hope that the analytic form of the strong coupling expansion of f{g) could be 
obtained from a perturbation theory for the matrix equation. 

Therefore, when formulated in terms of the Zmn, the problem becomes amenable to 
numerical study at all values of the coupling. We find that the numerical procedure 
converges rather rapidly, and truncate the series expansions of s{t) and of the kernel after 
the first 30 orders of Bessel functions. 
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Figure 1: Plot of the solutions of the integral equations: fm{g) for the ES kernel K^"^^ 
(upper curve, red), fo{g) for the kernel (middle curve, green), and f{g) for the 

BES kernel K^""^ + 2if (1 ower curve, blue). Notice the different asymptotic behaviors. 
The inset shows the three functions in the crossover region < g < 1. 

The function f{g) is the lowest curve plotted in Figure 1. For comparison, we also plot 
fm{g) which solves the integral equation with kernel K^™"^ [20], and foig) which solves 
the integral equation with kernel K^"^^ + K^^\ Clearly, these functions differ at strong 
coupling. The function f{g) is monotonic and reaches the asymptotic, linear form quite 
early, for g — 1. We can then study the asymptotic, large g form easily and compare it 
with the prediction from string theory. The best fit result (using the range 2 < g < 20) is 

, , , , , 0.0232 ±0.0001 , , 

f{g) = (4.000000 ± 0.000001)^ - (0.661907 ± 0.000002) + . . . (24) 

9 

The first two terms are in remarkable agreement with the string theory result (2), while 
the third term is a numerical prediction for the 1/g term in the strong coupling expansion, 
which may perhaps be checked one day against a two-loop string theory calculation. 

We do not need to restrict the numerical analysis to real values of g; complex values 
of g are of interest as well. In [19] it was argued that the dressing phase has singularities 
ai g ^ ±m/4, for n = 1,2,3, .. .. Also, their analysis of the small g series shows that 
there are square root branch points in f{g) at g = ±i/4. Perhaps, this is related to the 
cuts in the giant-magnon dispersion relations [33, 34, 35, 36, 27], for momenta close to 
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TT. Our numerical results indeed indicate branch points at ~ ±i/4, ±i/2 with exponent 
1/2. Beyond that we observe oscillations of both the real and imaginary parts of f{g) for 
nearly imaginary g. Further work is needed to elucidate the analytical structure of f{g)- 

3 Discussion 

A very satisfying result of this paper is that the BES integral equation yields a smooth 
universal function f{g) whose strong coupling expansion is in excellent numerical agree- 
ment with the spinning string predictions of [7, 9]. This provides a highly non-trivial 
confirmation of the AdS/CFT correspondence. 

The agreement of this strong coupling expansion was anticipated in [19] based on a 
similar agreement of the dressing phase. However, some concerns about this argument 
were raised in [22] based on the slow convergence of the numerical extrapolations. Luckily, 
our numerical methods employed in solving the integral equation converge rapidly and 
produce a smooth function that approaches the asymptotics (2). The cross-over region of 
f{g) where it changes from the perturbative to the linear behavior lies right around the 
radius of convergence, gc — 1/4, corresponding to gyu^ — For N — 3, this would 
correspond to ~ 0.25. 

The qualitative structure of the interpolating function f{g) is quite similar to that 
involved in the circular Wilson loop, where the conjectured exact result [31, 32] is 

The function (25) is analytic on the complex plane, with a series of branch cuts along the 
imaginary axis, and an essential singularity at infinity. The function f[g) is also expected 
to have an infinite number of branch cuts along the imaginary axis, and an essential 
singularity at infinity [19]. 

Let us compare this with the exact anomalous dimension of a single giant magnon of 
momentum p [33, 34, 35, 36, 27]: 

-1 + ^/i^ie^^hT^lj . (26) 

This function has a single branch cut along the imaginary axis, going from g — ^^^^-^ 

to infinity, and no essential singularity at infinity. Observables of the gauge theory are 
composites of giant magnons with various momenta p G (0, 27r). We thus expect the 
anomalous dimension of a generic unprotected operator to have a superposition of many 
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cuts along the imaginary axis. This should endow a generic multi-magnon state with an 
analytic structure similar to that of f{g)- 

We found numerically the presence, in f{g), of the first two branch cuts on the imaginary 
axis, starting ai g — n — 1,2. The first of them, which also occurs for the giant 
magnon with maximal momentum p = tt, agrees with the summation of the perturbative 
series [19]. The full structure of f{g) is expected to contain an infinite number of branch 
cuts accumulating at infinity, where an essential singularity is present. 

It is remarkable that the integral equation of [19] allows f{g), which is not a BPS 
quantity, to be solved for. Hopefully, this paves the way to finding other observables as 
functions of the couphng in the planar Af — A SYM theory. 
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